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Abstract 

The system of our interest is a dilute binary mixture, in which we consider that the species 
have different temperatures as an initial condition. To study their time evolution, we use the 
full version of the Boltzmann equation, under the hypothesis of partial local equilibrium for 
both species. Neither a diffusion force nor mass diffusion appears in the system. We also 
estimate the time in which the temperatures of the components reach the full local equilib- 
rium. In solving the Boltzmann equation, we imposed no assumptions on the collision term. 
We work out its solution by using the well known Chapman-Enskog method to first order 
in the gradients. The time in which the temperatures relax is obtained following Landau's 
original idea. The result is that the relaxation time for the temperatures is much smaller 
than the characteristic hydrodynamical times but greater than a collisional time. The main 
conclusion is that there is no need to study binary mixtures with different temperatures when 
hydrodynamical properties are sought. 
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1 Introduction 



Binary mixtures in which the temperatures are equal have been treated in literature for dilute 
gases [I] [SJ [3] , both in the case of inert components as well as those composed by charged particles 
thus leading to the study of plasmas [H [5J IE1 [Z] • These have been studied with kinetic theory when 
temperatures of ions and electrons are equal [SJ EU [TOJ QT] , and other authors. 

The question about the time relaxation for a binary mixture at different temperatures has been 
treated by Landau [TJl [13]. He used a Fokker-Planck kinetic equation. In this work we use the 
complete Boltzmann equation with no assumptions on the collisional term. We use the Chapman- 
Enskog method instead Grad's method because it is better defined in physical terms [14j[T5]. This 
allow us to give a microscopic foundation about the relaxation time. Such a calculation can be 
extended to plasmas but in discussion section we will give a few arguments that makes the exercise 
unnecessary. 

When we solve the Boltzmann equation with the Chapman-Enskog expansion to first order in 
gradients [TJ there is no diffusion force appearing in the linearized Boltzmann equation. This result 
deserves more discussion because it has generated misunderstandings in the literature [16] . 

This paper is divided as follows: Section 2 is devoted to some definitions and a review of kinetic 
theory as it relates to the problem. In section 3 we discuss the solution method and exhibit the 
linearized Boltzmann equation. In section 4 the time in which the temperatures reach the full local 
equilibrium is estimated. Section 5 is devoted to some concluding remarks. 

2 The kinetic theory scenario 

We consider a dilute binary mixture of particles with masses m a and mj,, The density is n = n a +ni,, 
so the mass density is p = p a + pb = m a n a + m^rtf,. Each species has a different local temperature, 
Tj where i = a, b. The Boltzmann equation in the absence of external forces for the evolution of 
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the distribution function for each species is 



in 



i,j=a 



where 



Afifj) = J - J [/(vD/K)-.f(v 4 )/(v,)] 



(2) 



xcr (viVj v-Vj) g lJ dv J dv' i dv' j 



where gij is the relative velocity between two particles when they collide. We remind the reader 
that there is a Boltzmann equation with the same structure for each species, coupled by J(fifj)- 
The cross section a satisfies the principle of microscopic reversibility guaranteeing the existence of 
inverse collisions. 

Next we introduce some useful definitions. In general the average of any dynamical variable ipi 
is given by 




(3) 



as usual we have the local particle densities 




(4) 



The local barycentric velocity is given by 




(5) 
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where the local velocity of each species is 

Ui(r,i) = — f fi(r,Vi,t)vidvi, (6) 
n t J 

implying that the chaotic velocity is defined as Ci = Vj u 4 . Recall that here we are taking u 4 as 
a local variable, because our hypothesis is that each species is in local equilibrium by itself. This 
is rather crucial in the calculation of the transport properties of the mixture since the mass flux 
Jj = it for i = a, b. The barycentric velocity u is not a candidate to be chosen as a local variable, 
because the temperatures of each species is different. 

The kinetic part of the stress tensor for the species i is defined as 

Ti = mi I fiCiCidvi. (7) 



Also, the heat flux for each species is given by: 

1 



J 9i = 2^( c » c i)' ( 8 ) 



and the mass flux is given by 



Ji=Pi(cj)=0 (9) 

meaning that there is no self-diffusion because the local velocity that we have taken is not the 
barycentric velocity of the mixture. Here Cj is the so called thermal or chaotic velocity. The 
average internal energy for the species i is, 

Pi£i = \pi{<%), (10) 
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which allows the introduction of the corresponding kinetic temperatures, namely, 

3 

= -UikTi. (11) 

Recall that the averages that represent macroscopic quantities are performed over Cj and not over 
Vi . That is because Vj = C; + U; includes the dynamical variable when the gas moves as a whole 
and in kinetic theory we desire to isolate the chaotic part of the velocity that generates thermo- 
dynamical fluxes [T71 [H] . 



To derive the conservation equations we take Eq. (1) plus the collision conserved quantities 
and integrate over the velocities c,. For the matter density we multiply by rrii to get, 

^ + V ■ (p^) = 0, (12) 

the expression for mass conservation. For the momentum balance, or Newton's second law for 
fluids, we multiply by mjVj, to get 

d 

— [piUi] + V ■ [n + p 4 u 4 Ui] = 0. (13) 

And for the balance for the internal energy density we multiply by ^rriicf 

dei f 1 

Pi-gj- + V • J gi + PiUi ■ Vet = / -m l f l c 2 i J{f l f j )dc i , (14) 

where i = a, b. If we sum over i making T a = Tf,, and substitute by u as a local variable we 
recover the mixture expressions under full local equilibrium. 
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3 Solution Method: Chapman-Enskog Expansion 

Since the Boltzmann equation for each species Eq. (1) has the same structure as the one in 
which the temperatures are not different, there is no need to undertake the details about its 
solution, which is given in Refs. [TJ [3J. The species are independent, as it can be seen from 
the balance equations. Thus, when we assume the validity of the functional hypothesis, namely 
fi (r, v|nj(r, t), Uj(r, t), ej(r, <)), we are explicitly indicating that the five independent state vari- 
ables for each species are rii, u», e^; i = a,b. Next fi may be expanded in a power series of Knud- 
sen's parameter e around the local equilibrium distribution function fjy , namely, the well-known 



Maxwellian distribution, 

ff 0) = n, (, 



\2^fj cxp r^r (5) 

These are the basis of the Chapman and Enskog method. To first order in gradients, we get that, 

fi=fP\l + <Pi), (16) 
so that from Eqs. JT]) and (J2j we have that, 

where the linearized collision kernels C(ipi) and C{tpi + <fj) are given by 

c (<Pi) = J J jfiidv' i - L dvidvi 1 a i igi t i 1 A(ipi), (18) 

where as in the one component case, the subscript 1 labels one of the two identical species. Also, 

C(<Pi +fj)= [ [ [ ./•"'/v;,/v, ( /v>,, // , ; A( 7 ;, + w), (19) 



G 



with 



A((fi) = $ + tp' n -tpi- (fa 



(20) 



and 



(21) 



Due to the functional hypothesis, the left hand side of Eq. (|T7|) can be expressed as 

ajf } _ d/| 0) dm £>/| 0) dm dfl 0) dTi 
dt dn t dt du, ' dt dT, dt ' 



(22) 



and the same for the other component. By using Eqs. (|15|) . (|12l) . (fT3|) and (fT4|) . as well as 
Pi = riikTi, the linear Boltzmann equation Eq. (|17|) turns out to be 



fT ci Ci : Vu 4 + --j Ci -^r L = [C{ Vi ) + Cfai + <pj)l (23) 

where c/cj = CjCi — (l/3)cfl. Here it is important to underline that a diffusive term which would 
be proportional to Vrij does not appear. This is a consequence of taking = n^fcTj instead 
p = p a + Pb = (n>a + n b) kT. Solution of Eq. (|2"3"j) is taken as the sum of the homogeneous plus 
inhomogeneous part [TH1 [2D] , 

= X • VlnT 4 +^ • VlnT, +^ : Vu; + B^j : Vu^ (24) 

1 2 
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Using now the subsidiary conditions jTTJ [22] 



AO) 



( \ 

rrii 



dci = 0, 



(25) 



in what follows, we shall omit the term Bi : Vu^ since by Curie's theorem [3T] does not couple 
with a thermal gradient. Therefore the solution to Eq. (|2"5]l is given by, 



<Pi = -tt^AVAT - -%AiVAT 



AT 



AT 



{ciAi +CjAj) VlnAT, 



(26) 



where AT = Tj — Tj. The still unknown functions A% (|cj|, rii, Ti) may be expanded in terms of the 
Soninc-Laguerre polynomials, 



oo 

Ai = 2_j a A^3/2 ( C i) 
p=l 



(27) 



where all properties of the coefficients a p A have been thoroughly discussed [3l |T0l E] and there is 
no need to undertake further unnecessary details. 



We now assume that the mixture of our interest is composed by electrons and heavy ions (pro- 
tons) with initial temperatures T® and T® respectively. Such an assumption allow us to consider 
that there is a huge difference between their masses. We consider that the ions are colder than 
the electrons because of the mass difference. Thus T° > T9, the superscript nought indicates their 
initial values. 



With the results here obtained we are now ready to face the problem. At some initial time t , 
T° > T® , there is a temperature gradient in the system which is as a whole in a non equilibrium 
state. Heat will flow through collisions from the mixture of electrons to that of the ions until after 
a certain time r elapses, T e = T p — T and the whole mixture is in a full local equilibrium state 
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with a certain temperature T. The question is, given and T® what is the magnitude of r? This 
will be discussed in the following section. 



4 Relaxation Time of the Temperatures 

Recall that the mixture is composed by electrons and protons whose the mass ratio is about 
1863, and rename the subscript as i — e for electrons and j = p for protons. The only way for 
the electrons to interchange energy with protons is when they collide. Then, due to the large 
difference of momentum and energy in a collision, small particles just change the sign of their 
velocities while large particles remain practically unaltered. This is the basis of Landau's idea, so 
that the basic mechanism for the thermalization are the electron-proton collisions. Whence, by 
using the definition of kinetic temperature, 

7^kT e = /im e Cg\ = ^- J m e c 2 e f e dc e , (28) 

we can write 

where AT = T e — Tj > 0. Further we have neglected the terms corresponding to the gradients 
in the left hand side of the Boltzmann equation, since when using the Chapman and Enskog ex- 
pansion contribute only to the second order. We also use the assumption of total ionization i.e. 
n e — n p = (l/2)n. 



Substitution of Eq. (JTJ in J2§J yields, 



|AT= ^JJJJ(m e c? e -m P 4){f>f>-fef P } (30) 
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Next, we take a collision model considering electromagnetic interactions Ref. [TT], so as shown in 
this reference, we have that, 



1 

n t 



C = J J CTdc' e d C ' v 



(31) 



(kT°) d 



hi 



1 + 16tt 2 



that has units [length] 3 /[time]. The initial temperature of the electrons is T® (notice that it has 
to be also substituted in fe°^). 

With the help of Eqs. ([26| and (|16l) . after substitution in Eq. (30) one obtains that, (see appendix) 



d_ 

dt 



AT = -h\VAT\ 



(32) 



with 



h=-™Cnm e 
k 




AT 



(33) 



where a\ is the coefficient that comes from the Sonine-Laguerre expansion and it is well-known 
[TO] to be of the order of the collisional time r ~ To derive last equation we had to use 

two approximations, the first one is that in a microscopic level, when one electron collides with 
a proton, the energy of both apparently remains unchanged because of the large mass difference, 
so c e ~ — c' e and c p ~ c' p . The second one is that we expanded in Taylor series up to first order, 
the logarithm In AT around the initial temperature difference AT . An outline of the calculations 
leading to Eq. (|33[) are given in the appendix. 



In order to integrate Eq. ([52)1 for simplicity we assume that the gradient VAT is in the x 
direction. The resulting equation is easily solved and the result is, 



AT(t) = AT°exp [-tit], 



(34) 
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note than in Eq. (I34[) the x dependence has disappeared due to the very particular boundary 
conditions. This is precisely what one expects, after a very short time the two temperatures 
must become the same at every space point in the mixture. The difference AT decays essentially 
as Landau's predictions. Numerically, by using for instance Tj? = 450K, T° = 200K and n ~ 
10 21 (l/m 3 ), the relaxation time is about 10~ 6 seconds. The relaxation time is not so sensitive to 
the initial difference AT , so if we insert appropriate temperatures (we mean, non relativistic) the 
order of magnitude does not change. This result shows explicitly that when we use kinetic theory 
to first order in the gradients, the relaxation time is small, and it is far from being relevant in the 
hydrodynamic regime. 

5 Conclusions 

We have calculated the relaxation time using a method which is entirely different as the one used 
by Landau [13] . based upon the Fokker-Plank equation. Here we used the complete information 
contained in the Boltzmann equation. In the Fokker-Plank approach the collision term in the 
Boltzmann equation is substituted with an "effective term" ; the problem in such equation is that 
it losses the microscopic information due to the collisions. In our method, as we have said, we 
imposed no assumptions on the collision term. In fact, the collisions are studied directly to find the 
relaxation time. This calculation allows us to justify form a microscopic point of view the hypoth- 
esis made in [TTj and many other works, that assume that in a binary dilute mixture or ionized 
plasma the temperatures are equal. We recall that this time is of the order of microseconds, so 
this process is very fast compared with those of hydrodynamic interest. To extend this calculation 
to a plasma, one needs to incorporate a term in Boltzmann's equation in which the Lorentz force 
is taken into account. Examining the solution in this case [HIE] when weak electromagnetic fields 
are considered, wont affect the order of magnitude of the relaxation time. 

It is also very important to discuss why the diffusive force docs not appear in these calculations 
and why no mutual diffusion is observed. When the temperatures are equal, a diffusive force 
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appears proportional to the gradient of the number of particles n, Ref. [35]. Now, since T e ^ T p , 
from our definition of chaotic velocity Cj = Vj + we have (cj) = 0, implying that the mass flux of 
each species is zero; we do not have mutual diffusion neither for the components nor the mixture. 
It is really important to highlight this point, the species are independent since they have totally 
independent maxwellian distributions. When the temperatures differ, an electron-proton collision 
produces a very small change in the energy of each one due to the large difference of the masses. In 
other words, in a collision between a very heavy particle and a very light one, the energy of them 
is almost unchanged. Only when the temperatures of the species become similar, we expect to see 
the diffusive effects. Thus, the diffusion will appear when the temperatures of the components are 
equal. 

A question could be raised about the applicability of this method to study relaxation times in 
the case of heavy Brownian particles suspended in a fluid. Although the question is meaningful it 
cannot be answered by the use of a Boltzmann type equation. In fact, the heavy Brownian particles 
undergo erratic motion due to the effect of many collisions between the particles of the fluid and 
them. That is why the process is better described by stochastic equations. If one would seek to 
find the relaxation time that takes to the Browninan particles reach the Maxwellian distribution, 
one have to solve such equations. For instance, a Langevin type or a Fokker-Planck type equation. 

Appendix 

The Chapmann-Enskog expansion to first order in the gradients is 

/ a = /i 0) +^, (35) 



where 



= -^AVAT - ^LA b VAT = - (c a A a + c b Ab) V In AT (36) 
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and AT = T a — T b . 

The expansion of the bilinear terms in /j appearing in equation in Eq. (30), with the use of 
Eq. (26) leads to the expression 



p ft _ f f _ AO)' AO)' _ AO) AO) 

JaJb Jajb — Ja J j J a J {, 



(37) 



= /i 0) / h (0) K + ^ - Va <Pb) 



AO) AO) 

J a Jfr 



- {c' a A a + c' b A b ) 

~ {c' a A a + c' b Ab) 

+ (c a A a + c b A b ) 
y + {c a A a + c b A b ) J 



V In AT. 



Using now the fact that protons are heavy particles (species b) and electrons light particles (species 
a) it follows that, in a collision c' a ~ — c a and c'h ~ c/,. Thus, the relative velocity g a b = |c a — c;,| ~ 
c a . Eq. (38) together with Eq. (30) leads to, 



d_ 

df 



AT = 



_4 



-C//K4-m fcC g)/i°Vr 



(38) 



x (4c a A a ) c a dc a dc b ■ V In AT. 



Direct integration in Eq. (39), using the Sonine-Laguerre expansion, the kinetic temperature 
|fcT b ° = ^rnb(cf) and the fact that rii = J f^dci leads to 



at k 



rpO 

' 5 1 5t 



a* |Vln AT | 



(39) 



which is Eq. (13"2"1) . 
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